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An investigation has been performed to determine which 
mathematical algorithms should be used in the calculation of atmospheric 
optical parameters using the Multlspectral Resource Sampler (MRS) sensor. 

A simulation of the MRS sensor was performed using a radiative-transfer 
model. The simulation provides the spectral radiance at the satellite 
sensor in terms of various atmospheric parameters, such as optical thickness, 
solar zenith angle, nadir view angle, relative azimuth angle, bi-directional 
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reflectance of the target, background albedo, and wavelength. 

Atmospheric correction algorithms were also developed for the 
determination of the total spectral optical thickness of the atmosphere 
for: 1} homogeneous (horizontal) hazy atmospheres, with diffuse targets; 

2) Inhomogeneous (horizontal) hazy atmospheres with diffuse targets; and, 
for, 3) homogeneous (horizontal) hazy atmospheres with non-diffuse targets. 
It was demonstrated, with numerous plots, that If the gonlometric (angular) 
properties of a neighboring pair of targets are known then the spectral 
optical thickness of the atnosphere can be detennined using the MRS sensor. 
Conversely, It was also shown that If one knows the optical properties of 
the atmosphere for a specified area then the bi-directional reflectance 
distribution function of the targets can be dv'termlned. 

As a result of this Investigation, It was concluded that the MRS 
sensor can be a valuable new device for the multlspectral remote sensing 
of Earth's surface. Specific recoonendatlons were made for more detailed 
Investigations of Inhomogeneous atmospheres and general bi-directional 
reflectance distribution functions. 
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PREFACE 


Tht Multi spectral Resource Sampler (MRS) "Proof-of-Concept" 

Study 1$ Intended to be a comprehensive analysis of the corrections that 
must be applied to MRS data to allow for atmospheric correction factors 
and the variability of bidirectional reflectance from the scene. 

This study was Initiated by Dr. Charles Schnetzler of NASA Goddard 
Space Flight Center, and was performed by ORI, Inc. Space Data and Systems 
Division with Mr. Charles W. Altken coordinating the efforts of ORI's 
consul tants . 

The complete study results are reported in five separate volumes 
which have the following titles and authors: 

DETERMINATION OF ATMOSPHERIC OPTICAL PARAMETERS USING 
THE MULTISPECTRAL RESOURCE SAMPLER, by Or. Robert E. Turner, 
Science Applications, Inc. 

ATMOSPHERIC CORRECTION USING AN ORBITAL POINTABLE IMAGING SYSTEM 
By Or. Philip N. Slater, University of Arizona 

A SIMUUTION ANALYSIS OF BIDIRECTIONAL REFLECTANCE PROPERTIES 
AND THEIR EFFECTS ON SCENE RADIANCE- IMPLICATIONS FOR THE 
MRS, By Or. James A. Smith, Colorado State University 

MRS LITERATURE SURVEY- OF BIDIRECTIONAL REFLECTANCE, By 
Or. James A. Smith, Colorado State Uhlversity 

MRS LITERATURE SURVEY OF ATMOSPHERIC CORRECTIONS 
By, Or. Philip N. Slater, University of Arizona 
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FOREWORD 


Th« NASA-designed Multlspectral Resource Sampler sensor 
is an entirely new device for the multispectral remote sensing 
of Earth's resources. It has many advantages ever existing 
sensor systems but the most important advantage from the 
atmospheric correction point of view is its pointability . 

The sensor can be set to receive radiation from a particular 
target or pair of targets for many nadir view angles. This 
fact is particularly advantageous because it allows an in- 
vestigator to determine the aerosol optical thickness of the 
atmosphere. 

Science Applications, Inc. was selected as a subcon- 
tractor to ORI to perform satellite simulations of the MRS 
satellite radiances and to develop algorithms which can be 
used to determine the optical thickness of the atmosphere. 

The optical thickness is the most important parameter for 
atmospheric correction of remote sensor data. 

Dr. Robert E. Turner was the principal investigator for 
this project. He performed the mathematical analysis and 
devised the computer program for the simulation studes. He 
would like to acknowledge the assistance of Mr. Joseph L. 
Manning who ran the con^uter program and helped in the plotting 
of graphs. 

In the course of this work informative discussions were 
held with Professor James Smith of Colorado State University. 

We also acknowledge the. direction provided by Hr. Charles 
W. Aitken of ORI, Inc. 
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IMTROOUCTIOM 


Zn th« last tan ysars an anornous data basa of information 
has baan eraatad on tha spactral radiation of tha Earth's 
surfaca. In aarly yaars« most of tha data wara collactad by 
using fflultispactral sansors aboard aircraft r wharaas in tha 
latar years most of tha data wara acquired with sansors on 
spacecraft such as Skylab and lAndsat. Of particular concern 
to many users of tha data is tha image degradation or tha 
probability of misclassification which arises from atmospheric 
affects . 

Earth's atmosphere contains a sami-parmanant suspension 
of particulates called an aerosol. Tha aerosol ma.iifests 
itself in tha form of light haze, industrial air pollution, 
heavy haze, fogs, and clouds. It is this aerosol component 
which is tha important factor in tha correction of remote 
sensor data for atmospheric af facts. 

In this report we shall be concerned with particular 
radiometric satellite tachniqu'^.s which can be used with the 
proposed MRS sensor to account for atmospheric effects in 
remote sensing. We describe a computer simulation of the 
Multispectral Resource Sampler (MRS) sensor for specific 
targets. In addition, we describe algorithms which can be 
used to determine the spectral optical thickness of the 
atmosphere provided the reflectances of "standard" targets 
are known. 
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REMOTE SENSING OF THE EARTH'S SURFACE 


Ssnsors on satsllitss in Earth orbit art usually dasi^nad 
to dataet alaetromagnatie radiation in a nunbar of spactral 
bands spaeifieally chosan so as to naxiaisa tha signal to 
background ratio for objacts of intarast. Tha ultraviolat 
ragion ganarally is not usad baeausa ataospharic gasas absorb 
most of tha radiation at thosa wavalangths. X.ikawisa« wa must 
taka spacial cara in tha inftarad ragion baeausa thara ara 
strongly absorbing gasas in tha infrarad. In addition, tha 
graatast amount of solar-raflactad anargy is in tha visibla 
part of tha spactrua. For thasa raasons tha naar ultraviolat, 
visibla and naar infrarad ragions eura chosan for radiomatric 
sansors. For tha analysis of ataospharic affects using tha 
MRS sensor tha atmospheric working group decided on tha 
spactral bands; 360-400 nm; 490-510 nm; 630-690 nm; and 
930-950 nm. 

In this section wa shall be concerned with tha description 
of tha natural radiation field in tha atmosphere and how it 
can be usad to account for and to determine tha optical prop- 
erties of tha atmosphere. 

2.1 RADIATIVE TRANSFER THEORT 

Radiation enters tha atmosphere and is scattered and 
absorbed by gases and particulates. If one has a description 
of tha optical properties of tha atmospheric gasas and par- 
ticulates in terms of altitude over a particular target area, 
axid a knowledge of tha terrain, than, in principle, it should 
be possible to calculate tha spactral radiation field. A 
complete description does not exist however, and even an 
approximate description of tha optical properties of tha 
atmosphere is difficult to obtain, at least for a larga 
portion of the Earth's surface. Nevertheless, if wa have 
soma information regarding tha reflectance properties of 
specific targets, than we can determine a quantity called the 
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spectral optical thicknats of tha ataosphara. Tha prinary 
raason that wa ara now ahla to accoapliah this task is that 
tha MRS sansor can "lock onto" a targat and obsarva that tar- 
gat as tha nadir viaw angla ehangas. This was not possibla 
with pravious satallita sansors. 

2.1.1 INTERACTION PROCESSES 

As radiation antara a seattaring-absorbing nadiua thara 
ara savaral proeassas or intaraction machanisas to considar 
in a nathaaatical dascription of tha radiation fiald. Thasa 
all illustratad for a voluna alanant in Figura 1. Within tha 
eontaxt of radionatry tha first procass is that of no intar- 
action (zaro gain and zaro loss) for radiation straaming into 
a sansor. Tha saeond procass is tha loss of radiation by 
absorption. Third# wa bava radiation which is seattarad out 
of t:ia fiald of viaw# rasulting in a loss. Fourth# thara is 
a gain of radiation which rasults from tha seattaring of 
radiation into tha fiald of viaw. Fifth# radiation can ba 
amittad by soureas in tha fiald of viaw. 

Tha loss maehanisms attanuata tha radiation batwaan 
sourca and sansor and this can ba axprassad as tha trans- 
mittanea along a path s# i.a.: 

I f* I 

T(X#s) « axp j -J ica#s')ds'j (1) 

Hhara ic(X#s‘') is tha voluaa axtinetion coafficiant at wava- 
langth \ and point s along tha path. 

Tha gain maehanisms of seattaring and amission giva risa 
to what is rafarrad to as tha path radianca# Lp(x#*#a)# a 
quantity which is a function of wavalangth X# position s# and 
diraction fi. 
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For tho pxurpostts of thif invostlgation wo do not nood 
to ba concornod with the details of the optical properties 
of the atmospheric constituents. We only need to understand 
the general definitions. An important optical quantity is 
the spectral optical thichness, defined as: 


T(X,e) 


* sec B 


I 


KiX.Z') d2 


( 2 ) 


where z' is the vertical distance from the Earth's surface, 
and 6 is the zenith angle. It should be noted that optical 
thickness is usually defined as Tq(X, 0), along a vertical 
path instead of a slant path. For this investigation how- 
ever, with its special emphasis on view angle we will use the 
definition of Equation 2. 

2.1.2 REMOTE SENSING EQUATION 

The spectral raidance is the amount of radiation flowing 
in a specific direction per unit area per unit time per unit 
solid angle per unit spectral interval. In the model used in 
this report we use units of milliwatts per square centimeter 
per steradian per micrometer of wavelength. Let the spectral 

A 

radiance at the sensor be L(X,0,Q) where 0 denotes the optical 
depth ? which is measured from the top of the atmosphere down- 

A 

wards. Let L(X,Tg,n) be the spectral radiance of the target at 
the target, i.e. at t ■ Tq. The radiance at the sensor can 
then be written as 


L(x,o,r.) 


Lq (a,Tq,.. }T(X,Tq,I! 


) + 


Lp (A,TQ,ii) 


( 3 ) 
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Eq. 3 can be called the remote sensing equation. If the path 

A ^ 

radiance L.(X«T,n) and the transmittance can be measured 

P 0 

or calculated then a measurement of the sensor radiance 

A ^ 

LU,0,fi) yields the target radiance LqU^Tq/^^}. It should 
be noted however, that the surface radiance is not independent 
of the atmosphere because it depends upon the radiance inci* 
dent upon the target and that is a function of the slcy radiance. 

An illustration of the MRS geometry is depicted in Figure 2 
for three nadir view angles and three azimuth angles. 

2.1.3 EQUATION OF TRANSFER 

If we consider the gain emd loss mechanisms illustrated 
in Figure 1, we can write an equation for the transfer of 
radiation. The time- independent, three dimensional radiative- 
transfer equation [1] is: 


Q • 7L(X,r,fl) » -<(X,r)L(X,r,a) 


+ j p(?,n . iV)L(X,r,.V)d;V 


+ a(X,r)B(X,r) +S(X,r,fl). (4) 

where <,a, and 3 are the volume extinction, volume absorption, 
and volume scattering coefficients respectively, and p is the 
scattering phase function. The left hand side of the equation 
defines the spatial rate of change in the radiance at point 

A 

r(x,y,z) in the medium for radiation in direction n. The 
first term on the right hand side defines the loss due to 
scattering and absorption; the second term describes the gain 
due to scattering into the field of view; the third term 
describes the gain due to thermal emission; and the last term 
describes some additional source. 
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FIGURE 2. VARIATION OF NADIR VIEW ANGLE FOR THE MRS FOR A SPECI 



The general three-dimensional problem is a very difficult 
one to solve. We shall consider the simpler one-dimensional 
problem in which there is a variation in the vertical direc- 
tion only. The corresponding one -dimensional equation [11 
(neglecting source S) is: 


M 


dL 

3t 


L(x,u»«) 


Uq(t) 

”■ 





b‘‘)L(T»u‘'»<^‘*)du 'd(p 



:t)) 


Uin (T); B(t) 


(5) 


where u * cos 9 and o>q(t) is called the single-scattering 
albedo. For the MRS bands we can neglect the last term in 
Equation 5 because thermal emission is negligible. 

2.1.4 BOUNDARY CONDITIONS 

The complete radiation field in the atmosphere is de- 
termined if one solves Eq. 5 subject to the following boundary 
conditions : 


L(O,u»0) “wEq (u + Ug)5(<> - <^g) 


( 6 ) 


and 


L(Tq,u» 4) ■ UgEg P(u»b/ - UQ/«g)e 




2, 1 

I I u"p(u,<0,u',P')Lg(TQ,u',?')du'dp' (7) 
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where Eq is the extraterrestrial solar irradiance, 6 q is 
the solar zenith angle and is the solar azimuth angle. 

The function is the general, bi-directional re- 
flectance of the surface and is the sky radi- 

ance. The radiance L(Tq,u>4>) is the same as the radiance 
Lo(^rTo,n) in Eq. 3 for a specific target. If one knows the 
bi-directional reflectance for the surface, then the surface 
radiance can be calculated using Eq. 7 provided the sky 
radiance is known. 

For a very clear sky with little haze the integral term 
in Eq. 7 is small compared to the direct solar reflectance 
term. On the other hand, for a hazy or overcast sky the 
first term is negligible. In another section we will use 
Eq. 7 to analyze the atmosphere using MRS correction algo- 
rithms . 

For a perfectly diffuse or Lambertian surface, Eq. 7 
becomes 


L("q,U/^) 




( 8 ) 


where J is the surface albedo which is independent of di- 
rection and E(tq) is the total (soleu: plus sky) irradiance. 

2.2 A RADIATION MODEL 

A number of computational methods have been developed 
over the last ten years which are used to calculate the 
spectral radiance in the atmosphere. In almost all cases 
there is usually an excessive amount of computer time for 
the computation because of the large number of scatterings to 
be represented. Chandrasekhar [1] was one of the first in- 
vestigators who found an exact solution to the radiation 
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transfer problem with multiple scattering in a molecular 
(Rayleigh) atmosphere bounded by perfectly diffuse 
(Lambertian) surfaces. Today/ many other computational 
techniques are used to calculate the radiation field. A 
brief description of most of the current methods is given in 
a State-of>the-Art report [2]. One of these methods, de- 
veloped by Turner 13,4,5,6,7] makes use of an approximation 
in the radiation field to arrive at a closed-form solution 
for the spectral radiance. The input and output quantities 
for this radiative- transfer model are illustrated in Figure 
3. Using the results of Coulson et al. [8] which were com- 
putations based upon Chandrasekhar's analysis, we compared 
the calculated radiances with those based upon the Turner 
radiative-transfer model. These comparisons are illustrated 
in Figure 4,5, and 6 for various sun angles and surface 
albedos. Although the atmosphere is a pure molecular one 
in which the horizontal visual range is 336 km, the agree- 
ment is very good. We have also compared the model with 
experimental data on sky radiance by Ivanov (9) for clear 
skies. These comparisons are depicted in Figures 7,8, and 9 
for different sun 2 uigles and again, the agreement is very 
good. 

It is interesting to look at a polar plot of the radiance 
as seen by a ground observer (sky radiance) and the total 
radiance (surface plus path) as seen by an observer in a 
satellite. Using the radiative-transfer model Turner [10] 
calculated and plotted these radiances as a function of zenith 
and nadir view angle for various surface albedos and atmo- 
spheric optical thicknesses. Examples of the radiances are 
illustrated in Figures 10 and 11 for Rayleigh scattering. 

The only difference between the two plots is the optical 
thickness. In Figure 11 the optical thickness is ten tLmes 
that in Figure 10 and the sky radiance is much greater. If 
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Input Quantities 


Wavelength 

Altitude 

Pressure 

Spectral Optical Thickness 
Single-Scattering Albedo 
Atmospheric Temperature 
Single-Scattering Phase Function 
Ground Temperature 
Ground Emissivity 
Target Reflectance 
Surface Albedo 
Solar Zenith Angle 
Zenith (nadir) View Angle 
Azimuthal View Angle 


Output Quantities 


Spectral 

Spectral 

Spectral 

Spectral 

Spectral 

Spectral 

Spectral 

Spectral 

Spectral 


Transmittance 
Global Irradiance 
Sky Irradiance 
Solar Irradiance 
Solar Radiance 
Sky Radiance 
Path Radiance 
Beam Radiance 
Surface Radiance 


Atmospheric Radiative Transfer Model 
Developed by r.e. Turner. 
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nCURE 5. OEPENDENCE OF SKY RADIANCE 
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38.9°. Vlsial range > 336 km; wavelength > 
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FIGURE 6 . DEPENDENCE OF SKY RADI- 
ANCE ON SCAN ANGLE (IN THE SOLAR 
PLANE); SOLAR ZENITH ANGLE • 0^. Vista] 
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nCURE 7. DEPENDENCE OF SKY RADIANCE 
ON SCAN ANCLE (IN THE SOLAR PLANE); 
SOLAR ZENITH ANCLE > eS.S**. Visual 
ranse • 269 km; wavaltngth • 0.404 um; surlaca 
albedo • 4S% 
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FIGURE 8. DEPENDENCE OP SKY RADUNCE 
ON SCAN ANCLE (PERPENDICULAR TO THE 
SOLAR PLANE); SOLAR ZENITH ANGLE • 
6S.S^. Visual raiun ■ 269 kra: wavelength « 
0.404 iim; surface albedo ■ 4S% 



nCURE 9. DEPENDENCE OF SKY RADIANCE 
ON AZIMUTH ANCLE. Visual range • 261 ka: 
wavelength • 0.447 urn; surface albedo • 16^; 
solar zenith angle > 65.7^. 


15 






Fit. 10 . nauam uw tky nOMiei at • Inaoios of *itw aa(k Id tha totar ptaao for a 
Kayldfh atSMiphcio and thnt nifleetaaeaa. d. iai3*S*;r0iiOLl;Xiiaj5«im:ibtday iiJiuw 
21; and D it 0.0, 0^, and 040. 


10 * * 0 * 10 


Fit . Total ladiiaot and iky ladianoa aa a Anedoo of liaw ao(la to tka tolar plana for a 
Raylaipk anoipkan and thna raflactanoob itSTl’ito is 14;X la OJS ids; tha day is Juna 
21; and a is 04. 043. and 040. 


6 




one adds an aerosol component to the atmosphere the scatter** 
Ing is predominantly in the forward direction. This effect 
is noticed in Figure 12 in which the total optical thickness 
is 0.10 and tha Rayleigh Optical thickness is 0.098. If 
the aerosol component is increased so that the total optical 
thickness is 0.2, characteristic of an atmosphere with a 
horizontal visual range of eUzout 2S0 km, then the radiance is 
that depicted in Figure 13. Here the strong forward scatter- 
ing is evident as well as the typical increase in the back- 
ward direction. If we increase the optical thickness to a 
value of 1, representative of a horizontal visual range of 
about 44 km, then the solar aureole is very pronounced as is 
illustrated in Figure 14. It should be noted that the scale 
has been changed tO accommodate the larger radiance value in 
the forward scattering direction. All of the illustrations 
are for Lambertian targets and backgrounds and for atmospheres 
which are uniform in the horizontal direction. In another 
section we will deal with non-Lambertian targets and in- 
homogeneous atmospheres. 

It is also of interest to observe the path radiance and 
the sky radiance as a function of ti.me (or sun angle) and 
the nadir or zenith view angle. These radiances are depicted 
in Figures 15 and 16 for a hazy atmosphere. The back-scatter- 
ing peak is clearly evident in the path radiance and the solar 
aureole is obvious in the sky radiance. The effect of surface 
albedo and visual range on the radiances is illustrated in 
Figxires 17 and 18. The multiple scattering of radiation is 
much greater for the shorter wavelengths and low visual 
ranges than for longer wavelengths and clear skies. 

In conclusion, the present radiative-transfer model seems 
to be in agreement with the more detailed, exact mathematical 
methods and also with experimental data. We will make use of 
calculations based upon this model for the MRS simulation. 
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ka; vuiol rangt ■ 9 )cb; lurfaet albtdo lor froon 
v^UUeo. 


F!?!* 






* North %• ^ 


nCURE 16. DEPENDENCE OF SKY RADIANCE ON 
TIME AND SCAN ANCLE. Soutbtastim Michicao, 

1 Sopttahor 1971. Altltodt • 0 ka; wavtiongth ■ 0.55 
um; Tisial rwge ■ 9 ka: tarlaet albodo lor grots 
vcgotation. 








FTCURE 17 DEPENDENCE OF PATH RADUNCE 
ON VISUAL RANGE AND WAVELpiGTH. Solar 
anflo ■ 30®; scan angle • 0°: azimuthal 
angle (in the plane of the sun);altitude •1km; 
surface albedo for green vegetation. 



FIGURE 13. DEPENDENCE OF SKY RADIANCE ON 
VISUAL RANGE .^ND WAVELENGTH. Solar zemth 
angle • 30®; scan angle • 0®; altltui^ ■ 0 km; azi- 
muthal angle • 0® (In the plane of the sun); surface 
albedo for green vegetation. 
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THE MULTISPECTBAL RESOURCE SAMPLER 


The Multispectral Resource Sampler (MRS) is an experi- 
mental pushbrooffl scan sensor being developed by NASA for Earth 
orbiting spacecraft flight in the mid-1980's. It's major 
advantage over existing systems is it's capability to point 
at and fix on a particular object on the Earth's 
surface as the spacecraft move, in an orbit. In this section 
we will use the radiative-transfer model to simulate the 
spectral radiances at the sensor for various angles, wave- 
lengths, atmospheres, and targets. 

3.1 MRS SIMULATION 

A technical publication by Schnetzler and Thompson [11] 
described the Multispectral Resource Sampler and its detailed 
spectral characteristics. The MRS sensor characteristics 
are listed in Table 1. The spectral bands to be used in the 
MRS depend upon the properties of the materials being 
investigated. These are illustrated in Figure 19. For the 
investigation of the atmospheric properties, one needs to 
know the atmospheric optical thickness in order to calculate 
the transmittance and path radiance. In this section we will 
use various optical thicknesses to calculate the spectral 
radiances at the sensor. The optical properties of the 
atmosphere, the angles to be considered, and the spectral 
reflectance of the vegetative background are given in Tables 
2 and 3. The nadir view emgles are those which correspond to 
equal increments in secant 9, i.e. 1.0, 1.1, 1.2, ... . 

The optical thicknesses are those which correspond to the 
simple Elterman [12] exponential haze model characterized 
by a surface horizontal visual range. Thus, we will consider 
V « 2 km (heavy haze) , V « 13 km (moderate haze) , and V « 40 km 
(light haze) . The single-scattering phase functions and the 
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OABLE 1 


MRS Sensor Characteristics 


Spectral Range: 
Spectral Bands: 


Spatial Resolution: 
Swath Width/Modes: 


Radioraetrtc Sensitivity: 
Data Rate: 

Pointabllity: 

Speed of Pointing: 


0.36 uin to l.Q uffi 

4 arrays, each with 2000 detectors 

5 selectable filters/array 
Bandwidths ^ 20 ran 
PolarizatlorT filters 


15 meters max 


15 ion at 15 m (4 bands) 

at 15 m (2 bands) 
at 15 m (4 bands, 502 sampling) 
30 Ion at 30 m (4 bands) 

Approximately 0.52 NEap (8 bit) 

15 mikfla 


2 axes 

♦ 40® across track 
+ 55® along track 

30®/sec across track 
SO/ser along track 
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TABLE 2. ANGLES FOR MRS SIMULATION 


1 Solar Zenith Angle 6^ « 30* 

60* 


Azimuth Angle 0 ■ 0*, 90®, 

180® 


Nadir View Angle 6*0®, 2i 

1.62®, 33.56' 

39.72®, 

44.42* 

•, 48.19®, 5] 

L.32®, 53.9 

56.25' 

', 58.24®, 6( 

• 

o 

o 

e 

• 


TABLE 3. OPTICAL PROPERTIES FOR MRS SIMOLATIOt 





A - Wavelength (uni) 

Tq - Optical Thic)cness with Surface Visual Range V * 

P « Albedo (Diffuse Reflectance) of Vegetative Background 


Deimendjian ' s Haze Model L was used for the aerosol size 
distribution with a pjurticulate refractive index of 
m» 1.5 - 0.01 i and a single-scattering albedo of 0.97. 
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single-scattering albedo used in the calculation of radiances 
are for an index of refraction m • 1.5 - 0.01 i, that is« for 
a slightly absorbing atmosphere. The background surface 
albedo is that for green vegetation. 

3.1.1 TARGET REFLECTANCES 

The general radiance at the surface in the upward 
direction is given by Eq. 7. If the bi-directional reflectance 
distribution function p were known then one could 

calculate the surface radiance. Rather than perform this 
tedious numerical analysis we will consider the surface 
radiance to have simple radiation patterns. We will use the 
following surface reflectances in our MRS simulation: 


Pl{6) 

• 0.25 

(9) 

PjO) 

• 0.25 cos 6 

(10) 

pjce) 

■ cos 6 sin ^6 

(11) 


The functions are represented graphically in Figure 20. The 
first function, p^(9) is a diffuse (Lambertian) surface and 
the other two represent hypothetical surfaces. All functions 
have azimuthal symmetry. 

3.1.2 MRS SENSOR RADIANCES 

Using the three reflectance functions above and the 
radiative- transfer model we calculated the sensor radiance as 
given by Eq. 3. In Figxire 21 we plot the radianc*e< at a 
wavelength of 0.380 ym vs. the secant of the nadir view angle 
for the Lambertian target. There is a large peak in the 
radiance in the anti-solar direction for the extremely hazy 
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(V a 2 km) atmosphere and a smaller peak for the light haze 
(V a 40 km) . It should be noted that the target reflectance 
for the diffuse target is one and not 0.25 as in Figure 20. 

The second target was used for the same situation and the 
resulting radiances are plotted in Figure 22. As expected, 
the radiance decreases more rapidly with Increasing angle 
because of the cosine surface distribution. The third re- 
flectance function was used and its corresponding radiances 
are illustrated in Figure 23. Here it should be noted that 
' the path radiance is for sec 6 ■ 1 because of the value of 
zero for the reflectance function when 8 • 0*. 

There should be a large change in radiance values for 
different wavelengths. This is illustrated in Figure 24 
for the diffuse target and in Figure 25 for the third target. 
As expected the backscattered radiation at the anti -solar 
angle is much greater for the short wavelengths than for the 
long wavelengths. For the large angles, however, the change 
in the transmittance is greater than the change in path 
radiance and the opposite relationship is found. 

All of these curves can change radically if the sun 
angle is different. This solar angle effect is illustrated 
in Figure 26, 27, and 28 for the three targets. As expected, 
the maximum is attained for the backscattering direction for 
the anti-solar angle, in this case when sec 6 > 2. 

Finally, we consider the case in which the orbit of the 
spacecraft is normal to the solar plane. As the sensor 
passes over the target the radiance distribution with angle 
should be symmetric for a horizontally homogeneous atmosphere. 
This effect is illustrated in Figures 29, 30, and 31. 
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FIGURE 26. RADIANCE AT 0.380 MM VS. THE SECANT OP THE NADIR VIEW 
ANGLE IN THE SOIAR PIANE FOR A WHITE DIFFUSE TARGET 
AND A VEGETATIVE BACKGROUND. SOLAR ZENITH ANGLE » 60 










(mn - js - jUis/MUJ)*! 
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FIGURE 31. RADIANCE AT 0.380 \m VS. THE SECANT OF THE NADIR 
VIEW ANGLE PERPENDICULAR TO THE SOLAR PLANE FOR 
A "COSINE-SINE-SQUARED" TARGET AND A VEGETATIVE 
BACKGROUND. SOLAR ZENITH ANGLE = 30®. 
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ATMOSPHERIC CORRECTION ALGORITHMS 


In this section we will consider several algorithms which 
can be used to determine the spectral optical thickness of 
the atmosphere. It is the atmospheric optical thickness which 
is of major importance in the calculation of the radiation 
field. This factor was clearly illustrated in the MRS simu- 
lation plots in the previous section. 

4.1 HOMOGENEOUS ATMOSPHERE - DIEFUSE TARGETS 

The most elementary case is that in which the atmosphere 
is uniform in the horizontal plane and for targets which are 
Lambertian. Because we are mainly concerned with the variation 
in the nadir view angle for the MRS, we will suppress many of 
the other parameters in the notation. The radiances at the 
sensor for two neighboring small targets are the following: 






^ u. 


P^^‘^0' 


( 12 ) 


and 








(13) 


If the atmosphere is uniform and the angular field of view 
is small enough so that there is a negligible effect of one 
target on another via atmospheric scattering then the trans- 
mi ttances and path radiances are equal. If we now form a 
difference function from Eq. 12 and 13 we get 
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D(Tq,9) 




(14) 


The transmittance in Eq. 14 is given by 


T(To,6) 


e”‘^0 


sec 6 


(15) 


so that we get the following equation for the logarithm of 
the difference function for two diffuse targets: 


£n D(Tq, 9) » £n Dq(Tq) - Tq sec 9 (16) 


The left-hand side of Eq. 16 is a measurable quantity for 
each spectral band of the MRS and, the nadir angle 9 is 
always known for each corresponding radiance. Thus, if one 
plots the logarithm of the difference function vs. the secant 
of the nadir view angle, a straight line results, the slope 
of which provides a value of the spectral optical thickness 
Tq. It is not important what the surface difference function 
Dq(Tq) is because that only determines the intercept. For 
this reason we will plot the logarithm of the spectral trans- 
mittance of a very hazy (V « 2 km) atmosphere vs. the secant 
of the nadir view angle in Figure 32. The slope is greater 
for the shorter wavelength curve because of the longer optical 
thickness at that wavelength. 

If we now consider one wavelength, say the one at 0.380 um, 
we can plot the logarithm of the transmittances vs. the secant 
of the nadir view angle for three atmospheres. This is illus- 
trated in Figure 33 for uniform atmospheres with horizontal 
visual ranges of 2 km, 13 km, and 40 km. 
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4.2 INHOMOGENEOUS ATMOSPHERE - DIFFUSE TARGETS 


Th« question arises, how does the atmospheric trans- 
mittance vary with the secant of the nadir view angle if 
the atmosphere is not uniform? Figure 34 illustrates the 
situation in which a satellite passes over a target and re- 
ceives radiation from various angles and therefore from 
various parts of a spatially variable atmosphere. In this 
case we have illustrated a situation in which there is an 
increase in the aerosol content in a horizontal plane but in 
which it is homogeneous in any vertical plane up to some con- 
stant height H. It is assumed that the molecular component 
is horizontally homogeneous everywhere. It should be noted 
that atmospheric refraction is negligible for the angle to 
be considered. 

We can now consider a simple model of a haze layer. 
Figure 3S depicts the geometry of the optical properties of 
the haze. The arbitrary origin is at x « 0 and the target 
to be considered is located at x. We will consider the 
visual range V(x) to be represented as the following: 


V(x) - V(0) 


I dV 1 d^ 2 

— — X + — X + 

II dx 21 dx^ 


(17) 


Here we shall limit our discussion to a "lineer" haze, i.e. 
one which is given by 


V(x) • V(0) + 


V(£) - V(0) 
£ 


X 


(18) 


where V(0) is the visual range at the origin and £ is the 
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scale length such that V(2-} is the visual range at l. 

The optical thickness along the path s for the point x 
at the surface is given by 


Tq{x,6) ■ j K(x', 


z')ds 


sec 6 


^ <(x',z')dz". 


(19) 


The total volume extinction coefficient <(x'‘,z'‘) is then 

■ so * 

V(x') 


( 20 ) 


where Sq. 20 is the usual relationship connecting extinction 
coefficient with visual range and it is to be understood that 
Eq. 20 holds for the wavelength X ■ 0.55 um. The k^(x',z') 
tern is the aerosol volune extinction coefficient and 
<q{2') i* the Rayleigh volume extinction coefficient. Using 
relation 18 in Eq. 19 we get 


*Note: In 50 ■ 3.912. This arises from the xise cf a 2% 

contrast for the human eye. 
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Tq(X/ 0) ■ £n 50 sac 9 


f 


dz' 


lo V(0) * x‘ 

X 


+ (Tg - Kg(0)H)sac 9 


whara 


T 


G 


I 


KQlz')dz' 


and <q( 0) is tha valua of <q(z) at z « 0. According to 
Figura 35 


x" * X + 2 ' tan 6 


and using this ralation in Eq. 21 givas us tha following 


Tg(X,9) 


X sac 9 In 50 


. — In 

V(x) - V(0) tan 9 


1 ^ 


H\V(x) - V(0) 


X V(X) 



<Q(0)Bj sac 9 


( 21 ) 


( 22 ) 


(23) 


tan 9 


(24) 
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Equation 20 describes the total optical thickness at an angle 
6 for a "lineeur" haze. It is interesting to see what this be- 
comes as 6 *>- 0. Define 


A . H [V(x) - V(6)1 tan 6 
X V(x) 


(25) 


then/ as A ->■ 0 Eq. 24 becomes 


Tq {X/9) 


H sec 9 in 50 
V(x) A 




+ (xg - <q(0) hJ sec e, 


(26) 


or 


Tq(x,0) 


H £n 50 
V(X) 




•Cq(O) E. 


(27) 


For other wavelengths we must consider a model for the haze. 
The model devised by Eltarzian [l2] will be used here for the 
spectral dependence. Thus, we have for the general, wave- 
length - dependent optical thickness of an inhomogeneous haze 
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tQ(^#X/9) 


X s«c 9 in 50 I(X) 


[V(x) - V(0)] tan 9 I(Xq) 


in 


1 + H fv(x) - V(0)1 
X V(x) 


T-(X) - <-(0,X-) H 

^ ^ ° KXq) 


sec 9 


(28) 


where X^ « Q.S5 un and I(X)/I(Xq] is given in Table 4. 


TABLE 4. SPECTRAL FUNCTION I(X)/I(Xq) 


X ( urn) 

0.380 

0.500 

0.660 

i 

0.940 

Function 

1.4536 

1.0876 

0.799 

0.5876 


The value of ■ 0.012 We will set the height 

of the hate layer to be 0.8621 kas, a value which is consistent 
with a horizontal visual range of 2 kn. We then establish this 
height for hazes of all visual rzmges but adjust the aerosol 
densities so as to provide the correct optical thickness. 

We can now determine the logarithm of the spectral trans- 
mittance as a function of the secant of the nadir view angle 
for various hazes, i.e. for different aerosol horizontal 
gradients. First, we consider the case where 


tan 9 
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V(0) « 40 km 
V(x) « 2 km 


Case 1 


and then 


V(0) • 40 km 
V(x) « 10 km 


Case 2. 


both with several values of x. The logarithm of the trans- 
mittance is given by 


£zi Tq i X / Xf 6 ) • — TQ(X/Xf9)« 


(29) 


It is assumed that the target is at the position x for which the 
visual range is V(x) . In Figure 36 we consider a situation of 
varying x for one of the MRS bands. The solid line represents 
a change in visual range of 38 km in a distance of 30 km« or 
a gradient of 1.267 km/km. The large dash curve is for an x of 
50 km and the small dash curve is for x « 100 km. The trans- 
mittance decreases quite rapidly as the view angle increases 
beyond 0"/ especially for the large gradient situation. As 
X «B the curves should become straight lines as in the 
homogeneous cases. It is interesting to note that the maximiun 
transmittance does not occur at 6 * 0” but rather on the high 
visual range side of the plot. The slope changes rapidly near 
6 » 0® in a somew.hat surprising manner. This can be easily 
seen if we plot the s^une set of curves on a highly magnified 
basis as in 
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Figure 37. It is easy to demonstrate that the slope of the 
curves at the inflection point (6 « 0”) is infinite if x is 
finite. The optical thickness is larger for the shorter 
wavelengths theui for the longer wavelengths. This effect is 
not surprising and is illustrated in Figure 38 for x « 30 km. 

In Figure 39 we illustrate the effect of a increased visual 
range at the point of observation ^V(x) > 10 kmj but with . 
higher gradients than in the lower visual range (v(x) « 2 km) 
case. Finally, we can understand the behavior of the trans- 
mittance in terms of the distance x for the positive angles 6 . 
This is illustrated in Figure 40 for sec 6 » 1.0, 1.6, 1.9, 
and 2.0. As the distance x decreases the gradient increases 
and the transmittance decreases because of the large optical 
thickness. 

From this analysis one can see that the spectral optical 
thickness of Earth's atmosphere can be determined if the trans- 
mittance can be measured as a function of the view angle. 

4.3 HOMOGENEOUS ATMOSPHERE - NON-DIFFUSE TARGETS 

Here we consider the more general situation in which 
the target is non-Lambertian. The basic "correction" 
algorithm is that of Eg. 14, i.e. 

0(Tq, 6) - DQ(TQ,6)T(TQ,e) (30) 


where the intrinsic difference function Dq(Tq, 9) is now de- 
pendent upon the view angle 6. If one makes the assumption 
that the sky radiance is constant over the hemisphere then 
Eg. 30 becomes 
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FIGURE 40. LOGARITHM OP Til’ TRANSMITTANCE vs. THE VISUAL RANGE SCALE 
DISTANCE X FOR FOUR SECANTS OF THE NADIR VIEW ANGLE. 
WAVELENGTH = 0.380 |im. V(0) = 40 km; V<X) = 2 km. 


D(Tq,0) 31 


“o^o* 




T(T, 



Uo''?o^ “ " ^0' 


+ L^(Tq)T(Tq,9) - 02^(U,P) 


(31) 


where ^ 3 ('^q) is the constant sky radiance and the p^^(u>'$) are 
the directional reflectances given by 


Pi^(u,4>) 


2it 1 



( U , P , u 


/P")du'dp'‘. 


(32) 


In order to illustrate the variation in D(T^,d) with sec 9 
we will only consider the direct solar contribution, that is, 
the first part of Eq. 31. Fig;:re 41 illustrates the variation 
in the di.fference function with the secant of the nadir view 
angle for two Lainbertian targets and three visual ranjas. The 
solar zenith angle is 60* but is of little importance e.xcept 
insofar as the absolute magnitude of the function is concerned. 
What Figure 41 indicates is just the variation of the trans- 
mittance along the line of sight as a function of the secant 
function. In Figure 42 we display the difference function for 
a diffuse surface and a sxirface with a cosine variation. From 
these two figxires one can easily see a large difference in the 
set of curves for the diffuse pair of targets and the diffuse- 
cosine target combination. In Figure 43 we illustrate a 
variation in the difference function for a combination of 
a diffuse target and one with a cos 9 sin ^ e variation. The 
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variability hara is aspacially pronouncad for tha high 
visibility casa. Finally > in Figura 44 wa prasant tha 
diffaranea function for a combination of tha tvro targats with 
angular variations. 

From this analysis wa hava shown that thara is anough 
variability in tha angular propartias of tha targats such 
thnt if ona knows tha gaomatric charactaristics of a targat 
than tha spactral transmittanca of tha atmosphara can ba 
datarminad from a maasuramant of tha diffaranea function* 
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CONCLUSIONS AND RECOMMENDATIONS 


In this investigation we have used a radiative-transfer 
model to simulate the MRS spectral radiances for several 
atmospheres characterized by various degrees of haze. We also 
performed the analysis for diffuse (Lambertian) targets and 
hypothetical targets with quite different goniometric properties. 
From the resulting radiance diagrams it is obvious that there 
can be large variations in the radiances depending upon the 
solar zenith angle > the nadir view angle, the relative azimuth 
angle, the optical thickness, and the target reflectances. 

In terms of its importance to the community of users of 
remotely sensed data the analysis of the variation of the 
difference function and transmittance is of greatest utility. 
Here we found a totally unexpected result, that is, the non- 
linear change in the transmittance as a function of view angle 
for inhomogeneous atmospheres and also the non-uniform change in 
the difference function for the hypothetical targets with 
various goniometric properties. If one knows the goniometric 
properties of selected targets, then a measurement of the 
difference function by the MRS yields the spectral trans- 
mittance of the atmosphere. Conversely, if one knows the 
spectral characteristics of the atmosphere a measurement of 
the difference function allows one to detenaxne the difference 
in the angular properties of the two targets. In addition, if 
one of the targets is known then the other is determined. 

Polarization has not been included in the present study. 

It is generally believed that polarization is not important 
because of the depolarization characteristics of the atmospheric 
particulates, expecially for very hazy atmospheres. 

From this preliminary analysis we have learned many things 
concerning the change in radiance for inhomogeneous atmospheres 
and non-Lambertian targets. Before the MRS sensor becomes 
"operational" however, we suggest that several additional 
studies be undertaken so that necessary design modifications 
can be made in the instrximent prior to launch. 
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First/ one should produce more diagrams of the simulated 
radiances for realistic targets with goniometric and radiometric 
properties which are characteristic of known natural materials 
and also for atmospheres with various amounts of absorptive 
aerosols. In this way one can simulate the MRS radiances for 
different classes of materials and atmospheric air masses. 

Thus, one might/ for example have radiances for black soil in 
an urban area of vegetation in a rural area. Certain areas 
may serve as calibrated test targets for atmospheric correction. 

Second/ the variation in transmittance with view angle 
for inhomogeneous atmospheres was surprising. It is now 
important that the haze gradient be specified for known 
environments characteristic of urban/ rural, coastal, etc. 
areas. Measurements should be made to determine the haze 
gradients. Also, one should model hazes more general than the 
simple "linear visual range" haze considered in this report. 

The derivative or slope of the transmittance plots could also 
provide information on the optical thic.kness for various 
points in a given area. Hence, the derivative of the trans- 
mittance function should be plotted. 

Third, we have seen how information can be obtained on 
the goniometric properties of surfaces if we know something 
about the atmosphere. One should definitely use realistic 
bi-directional reflectance distribution functions with in- 
homogeneous atmospheres to calculate the difference function 
for the two target radiances for the MRS spectral bands in 
order to determine the sensitivity of the MRS radiances for 
various targets. 

If these tasks are undertaken then a definite assessment 
can be made of the potential uses of the MRS for atmospheric 
correction of remote sensor data. In addition, a whole new 
remote sensing technique may result for the classification of 
Earth resources based upon their goniometric properties as 
well as their spectral characteristics. 
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